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This  report  consists  cf  three  volumes  snd  documents  the  current  version  in 
the  C81  family  of  rotcrcraft  flight  simulation  programs  developed  by  Bell 
Helicopter  Company.  This  current  version  of  the  digital  computer  program  is 
referred  to  as  AGAJ7'\ - - 

The  new,  revised,  or  alternate  mathematical  models  incorporated  into  the 
program  dux 'eg  the  current  contract  are  as  follows: 
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Block  20.  Continued 

(1)  Fuselage  aerodynamic  forcea  and  momenta  (reviaed) 

(2)  Aerodynamic  aurfacea  (reviaed  with  two  aurfacea  added) 

(3)  External  atorea/aarodynamic  brakea  (new)  ^ 

(4)  Rotor  blade  airfoil  aectlon  dlatributlon  (new) 

(5)  Rotor-induced  velocity  dlatributlon  (alternate) 

(6)  Rotor  unateady  aerodynamica  (alternate) 

(7)  Rotor  wake  effect  at  aerodynamic  aurfacea  (alternate) 

(8)  Method  for  numerically  integrating  rotorcraft  equatlona  of 
motion  (alternate) 

Thla  volume,  the  Programmer' a  Manual,  containa  the  information  ntceaaary 
to  aet  up  and  aupport  the  computer  program.  Specifically,  it  includes  cross- 
ref  erencea  of  FORTRAN  COMMON  BLOCK  variables,  a  catalog  of  subroutines,  and 
a  discussion  of  programming  considerations.  The  listings  and  related  soft- 
were  for  the  computer  program  documented  in  thla  report  are  unpublished  data 
iriilch  are  on  file  at  the  Euatis  Directorate,  U.  S.  Army  Air  Mobility  Research 
and  Development  Laboratory  (USAAMRDL),  Fort  Euatia,  Virginia.  Volume  1,  The 
Engineer's  Manual,  documents  the  background  and  development  of  the  current 
version  of  the  program.  Volume  II,  the  User's  Manual,  contains  the  detailed 
information  necessary  for  setting  up  an  input  data  deck  and  interpreting  the 
computed  data. 
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This  report  and  its  accompanying  computer  program  were  developed  under 
Contract  DAAJ02 -72-C-0098  awarded  in  June  1972  by  the  Eustis  Directorate 
oi  the  U,  S.  Army  Air  Mobility  Research  and  Development  Laboratory 
(USAAMRDL).  In  addition  to  the  work  performed  under  this  contract,  the 
report  and  computer  program  include  the  documentation  and  program  features 
developed  under  USAAMRDL  Contracts  DAAJ02-70-C-0063  and  DAAJ02-73-C-0086. 
The  contractor  and  USAAMRDL  have  agreed  that  the  computer  program  docu¬ 
mented  herein  is  the  new  master  version  of  the  prograrfi.  Hence,  this 
report  supercedes  all  previous  versions  of  the  C81  program  and  documenta¬ 
tion. 

Technical  program  direction  was  provided  by  Mr.  E.  E.  Austin  of  USAAMRDL. 
Principal  Bell  Helicopter  personnel  associated  with  the  current  contract 
were  Messrs.  B.  L.  Blankenship,  J.  M.  Davis,  and  P.  Y.  Hsleh,  and  Dr.  B.  T. 
Uaak.  In  addition,  Dr.  R.  L.  Bennett  and  Mr.  B.  J.  Bird  assisted  in 
coordinating  the  work  and  documentation  prepared  under  the  two  previous 
contracts  noted  above  with  that  prepared  under  this  contract. 
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This  manusl  documents  the  rotorcreft  flight  simuletion  program,  designated 
AGAJ73,  and  its  post  processor  for  data  reduction,  designated  GDAJ07 .  To 
the  user,  this  system  appears  as  a  single  program;  to  the  programmer,  the 
two  programs  are  very  different.  This  documentation  is  fcr  the  programs 
as  they  were  written  for,  and  are  being  used  on,  an  IBM  System/360  Model 
65  Computer  at  Bell  Helicopter  Company. 


The  Information  in  this  volume  is  of  two  types.  Section  2  contains  the 
information  necessary  to  get  the  programs  operational  on. a  computer  com¬ 
patible  with  the  Mi  tallatlon  at  Ball  Helicopter  Company.  If  the  pro¬ 
grams  are  to  be  modified  in  any  way,  the  programmer  will  need  the  infor¬ 
mation  in  Sections  3  and  4  of  this  volume. 


I 


F 

0 


2.  OPERATING  ENVIRONMENT  AND  PROCEDURES 

The  System/360  environment  under  which  this  program  system  is  maintained 
is  Operating  System/360  Option  IV  (MVT).  Input  on  the  system  reader  is 
controlled  by  Houston  Automatic  Spooling  Program  (HASP)  II ,  as  is  system 
output  destined  for  on-line  printer  or  card  punch.  As  delivered  under 
this  contract,  the  program  requires  approximately  a  1000K  region  of  memory 
(IK  ■  1024  bytes,  1000K  «  1  megabyte).  Methods  for  reducing  this  large 
memory  requirement  are  discussed  below  and  in  Section  .4.4.  The  Basic 
Sequential  Access  Method  routines  are  retained  in  a  resident  LINKPACK 
area.  Scratch  data  sets  are  directed  to  a  Telex  5312  Direct  Access 
Storage  Facility.  Tape  data  sets  are  recorded  on  Telex  6420  Tape  Drives. 
The  recording  format  for  the  CALCOMP  PLOTTAPE  is  9  track,  800  bpi,  NRZI. 

All  other  tape  data  sets  are  recorded  as  9  track,  1600  fci  phase  encoded. 
To  run  the  program  under  0S/360,  at  le'?t  262,144  (256K)  bytes  of  main 
memory  are  required  for  operation  with  Option  I  (POP)  or  Option  II  (MFT) 
systems.  At  lesst  524,288  (S12K)  bytes  of  main  memory  are  required  for 
Option  IV  (MVT). 

The  program  has  been  maintained  on  the  IBM  System/360  FORTRAN  H  Compiler 
in  USA  FORTRAN.  Compiler  options  used  are  SOURCE,  EBCDIC,  NOLIST,  NODECK, 
LOAD,  MAP,  NOEDIT,  ID,  XREF,  and  OPT-2.  Since  the  FORTRAN  H  Compiler 
performs  essential  optimization  functions,  compilation  of  this  program 
under  any  other  compiler  or  using  optimization  level  less  than  "2"  on  the 
H  Compiler  will  result  in  degraded  performance  in  speed.  Some  FORTRAN 
language  features  peculiar  to  IBM  FORTRAN  definitions  have  been  noted 
explicitly  in  individual  routines. 

AGAJ73  can  be  link-edited  in  several  different  ways.  For  example,  with 
the  OVERLAY  structure  shown  in  Table  2-1, the  entire  program  can  be  loaded 
into  either  main  or  extended  memory)  with  the  HIARCHY  support  structure 
shown  in  Table  2-2, part  of  the  program  can  be  loaded  into  main  memory  and 
the  rest  into  extended  memory)  or  If  neither  OVERLAY  nor  HIARCHY  structure 
is  used,  the  entire  program  can  be  loaded  into  either  main  or  extended 
memories.  The  best  way  to  link-edit  this  program  depends  upon  the 
facilities  available  at  the  local  installation. 
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Although  this  program  is  quite  Urge,  it  can  be  run  in  less  than  300K  if 
the  OVERLAY  structure  shown  in  Table  2-1  is  used  and  the  size  of  several 
large  arrays  is  reduced  according  to  the  procedure  outlined  in  Section 
4.4.  If  the  region  size  available  is  greater  than  300K,  the  OVERLAY 
shown  in  Table  2-1  can  be  relaxed  to  reduce  the  overhead  time  associated 
with  OVERLAY.  By  using  the  HIARCHY  support  structure  shown  in  Table  2-2 
and  also  following  the  procedures  in  Section  4.4,  the  program  can  run  in 
a  206K  region  of  main  memory  and  a  450K  region  of  extended  memory.  It 
would  need  approximately  another  300K  if  the  reduction  of  array  sizes 
mentioned  above  is  not  done.  Since  the  buffers  are  in  extended  memory, 
the  amount  of  extended  memory  used  depends  upon  the  units  opened. 
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The  input  data  for  the  linkage  editor  for  GDAJ07  are  ahown  in  Table  2*3. 
This  OVERLAY  atructure  reaulta  in  a  program  whose  longeat  segment  is  142K 
bytea.  After  a  buffer  allocation,  the  size  of  the  region  needed  ia  184K 
by tea. 

The  Job  Control  Linkage  (JCL)  uaed  to  run  a  typical  aet  of  data  ia  shown 
in  Figure  2-1.  If  other  than  HIARCHY  structure  ia  uaed,  the  REGION  card 
in  the  first  atep  ahould  be  changed  accordingly  from  that  ahown  in  Figure 
2-4.  Tablea  2-4  and  2-5  give  the  unit  type,  what  it  ia  uaed  for,  and 
which  aubroutinea  refer  to  it  for  AGAJ73  and  GDAJ07 ,  respectively. 

The  JCL  has  been  prepared  such  that  the  Job  is  in  two  steps:  C81  and 
C81FL0T.  The  C81  step  executes  AGAJ73,  and  the  C81PL0T  atep  executes 
GDAJ07  even  if  the  first  step  abends.  This  was  done  ao  that  if  a  maneu¬ 
ver  abends  on  the  first  step,  the  user  will  still  get  plots  of  what  was 
run. 


TABLE  2-1.  LINKAGE  EDITOR  CARDS  FOR  OVERLAY  OF  AGAJ73 


OVERLAY  ALPHA  INITIALIZATION  SEGMENT 

INSERT  PDSREOf REOIO»REORWKf  REOSMK »  START  , WKTABN, IHCNAMEL 
OVERLAY  BETA 

INSERT  FPRCHK»JSTKED,LGCINT ,L I ZE ,NPUTOT ,RBDY,REAOIN 
INSERT  REOATBf RECBMSfRCOCL »  XST INT  ,  ZERO 
OVERLAY  BETA 

INSERT  FUSINT,  INBLD,  INRl),  INRTR»  I  NSC  AS  «  J  FBG I N, MODAL  •  PYL  I  NT 
INSFRT  RTTNIT*STBZIN|W9M0AL*XC0NIN 
OVERLAY  BETA 

INSERT  KNEM,T  ABF  IX t T ABOUT , TURN,YR INIT,YSINIT 
OVERLAY  ALPHA  HARMONIC  ANALYSIS 
INSERT  HARM.LOADT 

OVERLAY  ALPHA  STABILITY  ANALYSIS  SEGMENT 
INSERT  ALLMAT, ALSTAB, INVERS .IOMAT, MORORSjMODFSt PUNCH, HRMS 
INSERT  IHCCLABS , IHCCI A S, I HCCL  SOT , IHCCS ATS , I HCCSAS , I HCFOPT 
INSERT  IHCLSQRT 

OVERLAY  ALPHA  GENERAL  PURPOSE  SEGMENT 

INSERT  ANAL, ANOOIT, AZMI  NT,  AZMUTH,CDCL,CLCD,CMCALC, DIFFER 
INSERT  OOTX, FOCUS, FORCMC, FUSFNM, HRESP, INTFRQ, ITROT  ,MRAL 
INSERT  NSTEOtPAOBON,  RADIAL,  RAOOUT,RGlJST ,  ROT  AN,  RT  WAKE,  SOLVE 
INSERT  STBWAK,STBZFM, SWSR AT , TABINT r TrmANU,UNSTED, W1 NG , WRFM 
INSERT  WRTMNV»W$HDUF,X STORE , IHCSASCN 
OVERLAY  BETA 

INSERT  AJACOR, JACOBI, DAMPER, WRVP 

OVERLAY  GAMMA  STABILITY  ANALYSIS  SEGMENT 

INSERT  WRINST 

OVFRLAY  DELTA 

INSFRT  INSTAB 

OVFRLAY  DELTA 

INSFRT  STAB,WRSTA'.« 

OVERLAY  GAMMA  TPJn  SEGMENT 
INSERT  TRIM.TVT 
OVERLAY  DELTA 
INSERT  RPTPG 
OVERLAY  DELTA  ' 

INSERT  TRMi 
OVFRLAY  OFLTA 
INSERT  WRTRIM 
OVERI  AY  DELTA 
INSFRT  ITR I M*  POZ1 
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TABLE  2-1.  Concluded. 


OVERLAY  BETA  MANEUVER  SEGMENT 

INSERT  HER  I V«HAMM$  »HPCG»MANU*MANUV » OUAN  »RNGKTA»  SAVfc IC»WAG 

OVERLAY  GAMMA 

INSERT  SCASIT 

OVERLAY  GAMMA 

INSERT  VAR  I 

OVERLAY  DELTA 

INSERT  VCNTRl 

OVERLAY  DELTA 

INSERT  SUPFMP 

OVERLAY  DELTA 

INSERT  MTLT 

OVERLAY  OELTA 

INSERT  GUST 

OVERLAY  DEI  T A 

INSERT  FXTORS 

OVFRLAY  OFLTA 

INSERT  CNTM 

OVERLAY  OELTA 

INSERT  MOMB 

OVERLAY  ONE(REGICN) 

INSERT  IN1T  »RFSTRT  »SIVAR» TIMLP»T IVAR»WRMANU 

OVERLAY  ONE 

INSERT  CONS  TR»  STBO 

ENTPY  main 
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TABLE  2-2.  LINKAGE  EDITOR  DATa  CARDS  FOR  HIARCHY  SUPPORT  OF  AGAJ73 


HIARCHY 

1 

IMCCOMH2 

HIARCHY 

It  INSTAR 

HIARCHY 

SCASIT 

HIARCHY 

I 

IMCECONH 

HIARCHY 

1« INTERS 

HIARCHY 

SIVAR 

HIARCHY 

1 

IHCEPIOS 

HIARCHY 

ItlONAT 

HIARCHY 

STANAN 

HIARCHY 

1 

IHCEFNTH 

HIARCHY 

It JACOBI 

HIARCHY 

STARAN 

HIARCHY 

1 

1MCCRRM 

HIARCHY 

It JFBGIN 

Diarchy 

START 

HIARCHY 

1 

IHCCTRCH 

HIARCHY 

ItJSTREO 

HIARCHY 

STBMAK 

HIARCHY 

1 

IHCFCVTH 

HIARCHY 

ItLGCINT 

HIARCHY 

STBZIN 

HIARCHY 

l 

INC  FORT 

HIARCHY 

ItLIZE 

HIARCHY 

SUFERF 

HIARCHY 

l 

IHCNANEL 

HIARCHY 

ItlOAOT 

HIARCHY 

SMSRAT 

HIARCHY 

1 

IHCUATBl 

HIARCHY 

It MAIN 

HIARCHY 

TABFIX 

HIARCHY 

1 

IHCUQFT 

HIARCHY 

ItNANU 

HIARCHY 

TABOUT 

HIARCHY 

1 

ABDUMF 

HIARCHY 

ItNAHUY 

HIARCHY 

TILT 

HIARCHY 

I 

A  JACOB 

HIARCHY 

ItNORORS 

HIARCHY 

TINLF 

HIARCHY 

I 

ALSTAB 

HIARCHY 

ltNHEN 

HIARCHY 

TIVAR 

HIARCHY 

1 

COZARH 

HIARCHY 

ItHOOAL 

HIARCHY 

TOFLOT 

HIARCHY 

t 

CHTM 

HIARCHY 

ItNOORS 

HIARCHY 

TRIH 

HIARCHY 

1 

COHSTB 

HIARCHY 

ItMOMB 

HIARCHY 

TRMAMU 

HIARCHY 

1 

CONTRH 

KIARCHV 

ItHTLT 

HIARCHY 

TRHl 

HIARCHY 

1 

OATC 

HIARCHY 

ItHOFS 

HIARCHY 

TURN 

HIARCHY 

1 

OANFER 

HIARCHY 

ItHFUTOT 

HIARCHY 

VARI 

HIARCHY 

l 

OERf  V 

HIARCHY 

ItFOSRIO 

HIARCHY 

VCNTRL 

HIARCHY 

& 

ERRCHR 

HIARCHY 

ItFOZI 

HIARCHY 

NAG 

HIARCHY 

1 

EXTORS 

HIARCHY 

ItFUNCH 

HIARCHY 

MRTABN 

HIARCHY 

1 

FOCUS 

HIARCHY 

ItFYlIHT 

HIARCHY 

NRFN 

HIARCHY 

1 

FORMA 

HIARCHY 

ItOUAN 

HIARCHY 

MRINST 

HIARCHY 

1 

FORMA l 

HIARCHY 

ItRAOOUT 

HIARCHY 

MRNANU 

HIARCHY 

1 

FORYYY 

HIARCHY 

ItRBOY 

HIARCHY 

MRNDAL  * 

HIARCHY 

1 

FORVVI 

HIARCHY 

ItREAO'N 

HIARCHY 

URNS 

HIARCHY 

I 

FOSMA 

HIARCHY 

1, REOAT  A 

HIARCHY 

MROTI 

HIARCHY 

FOSMAI 

HIARCHY 

ItRFOBHS 

HIARCHY 

MR STAB 

HIARCHY 

1 

FUSIHT 

HIARCHY 

ItRCOCL 

HIARCHY 

MRTNNV 

HIARCHY 

1 

GUST 

HIARCHY 

ItREOlO 

HIARCHY 

MR  TRIM 

HIARCHY 

l 

HAIMS 

HIARCHY 

ItREORMK 

HIARCHY 

MRMF 

HIARCHY 

1 

HFCO 

HIARCHY 

ItREOSMK 

HIARCHY 

ICON IN 

HIARCHY 

1 

INBLO 

HIARCHY 

ItRESTRT 

HIARCHY 

XSTINT 

HIARCHY 

l 

INI? 

HIARCHY 

ItRNGKTA 

HIARCHY 

VRINIT 

HIARCHY 

1 

INRO 

HIARCHY 

ItRFTFG 

HIARCHY 

YSINIT 

HIARCHY 

1 

INRTR 

•  HIARCHY 

ItRTINIT 

HIARCHY 

ZERO 

HIARCHY 

1 

INSCAS 

HIARCHY 

ItSAVEIC 

HIARCHY 

ZLLCAL 

HIARCHY 

1 

INSTAB 

HIARCHY 

ItSAYOLO 

ENTRY  MAIN 

TABLE  2-3.  LINKAGE  EDITOR  DATA  CARDS  FOR  OVERLAY  OP  GDAJ07 

IE WO OOO 

OVERLAY  ALPHA 

IFWOOOO 

INSERT  CURVET 

IEWOOOO 

INSERT  IHCLATN2 

IFWOOOO 

INSERT  1HCLS0RT 

IEWOOOO 

INSERT  C81L 

IEWOOOO 

INSERT  TIMPTS 

IEWOOOO 

OVFRLAY  ALPHA 

IEWOOOO 

INSERT  CALCON  * 

IEWOOOO 

INSERT  IHCFRXP I 

IEWOOOO 

INSERT  LINE 

IEWOOOO 

INSERT  NEXTTIME 

IFWOOOO 

INSERT  NUMBER 

IFWOOOO 

INSERT  SYMBOL 

IFWOOOO 

OVERLAY  BETA 

IEWOOOO 

INSERT  FSFT 

IEWOOOO 

INSERT  HARM 

IEWOOOO 

INSERT  IHCSSOPT 

IEWOOOO 

INSERT  PLOTFR 

IEWOOOO 

INSERT  AXIS# 

IEWOOOO 

INSERT  IHCFRXPR 

IEWOOOO 

INSERT  IHCLSCN 

1  IEWOOOO 

INSERT  IHCSEXP 

IEWOOOO 

INSERT  IHCSLOG 

IEWOOOO 

INSERT  SCALE# 

IEWOOOO 

OVFRLAY  BETA 

IEWOOOO 

INSERT  SCALIT 

IEWOOOO 

!NSCRT  CALCB1 

I FWOOOO 

INSERT  PPLOT 

IEWOOOO 

INSERT  SCLFIX 

IEWOOOO 

INSERT  INPLOT 

IFWOOOO 

• 

ENTRY  MAIN 

\ 

\ 

\ 

\ 

l-  ,  _ 1  1  . -  .  ,  . =a 
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TABLE  2-4.  INPUT/OUTPUT  UNITS  USED  IN  AGAJ73 

Unit 

No. 

Typo 

Used  for 

Used  by  Subroutine 

1 

direct 

access 

Permanent  data  storage 
in  partitioned  data 
set  (PDS) 

JSTRED,  REDATB,  REDBMS, 

REDCL,  REDID,  REDRWK, 

REDSVK  . 

2 

tope 

New  restart  tape 

RESTRT 

3 

direct 

access 

Utility  storage  of  time- 
variant  trim  flapping 
history  and  storing 
maneuver  time  history 
to  pass  to  GDAJ07 

INIT,  LOADT,  MAIN,  MANU, 
MANUV,  RESTRT,  TVT  _ 

4 

tape 

Old  restart  tape 

RESTRT 

5 

card 

reader 

Input  data 

MAIN 

6 

printer 

Printed  output 

ALSTAB,  AZMUTH,  CDCL,  CLCD, 
DAMPER,  ERRCHK,  EXTORS, 

FUSINT,  HAMMS,  HPCG,  INRO, 
INSTAB,  IGMAT,  ITRIM,  ITROT, 
JFBCIN,  LIZE,  LOADT,  MAIN, 
MBAL,  MNEM,  NPUTOT,  RADOUT, 
READ IN,  REDID,  REDRWK,  RPTPG, 
SIVAR,  STAB,  TABOUT,  TIVAR, 
TRIM,  TRM1,  TURN,  VIND,  WAG, 
WRFM,  WRINST,  WRMANU,  WRMDAL, 
WRMS,  WROT1,  WRSTAB,  WRTMtV, 
VRTRIM,  WRVP,  XOONIN,  YRINIT, 
YSXNIT 

7 

card 

punch 

Punched  output 

PUNCH 

• 

10 

direct  . 

access 

Utility  storage  of 

AQAJ73  input  data 

JSTRED,  MAIN,  READIN,  REDATB, 
REDBMS,  REDCL,  REDID,  REDRWK, 
REDSVK  _ 

11 

direct 

access 

Passing  input  data 
to  GDAJ07 

MAIN,  READIN 

v. 
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TABLE  2-4.  Continued. 

Unit 

Mo. 

Typ# 

Used  for 

Used  by  Subroutine 

14 

direct 

eccees 

Stored  end  retrieved 
Inltlsl  conditions 
when  Henning's  method 
Is  used 

SAVEIC 

TABLE  2-5.  INPUT/OUTPUT  UNITS  USED  IN  GDAJ07 

Unit 

No. 

Type 

Used  for 

Used  by  Subroutines 

3 

direct 

eccess 

Maneuver  time  history 
from  AGAJ73  or  Tape  8 

CURVET,  C81L,  PSPT,  MAIN, 
SCALIT 

j  6 

printer 

Printed  output 

CALCS 1,  CURVET,  C18L,  PSPT, 
PPLOT,  WROT1 

8 

tspe 

Old  time  history  tape 

C81L 

9 

tape 

New  time  history  tape 

C81L 

10 

direct 

eccess 

Input  date  from 

AGAJ73 

CURVET,  C81L,  PSPT,  MAIN, 
SCALIT 

PLOT- 

TAPE 

tspe 

Plot  maneuver  time 
history  In  GDAJ07 

PLOTER 

2-8 


mis* 


//IfVIJtfO  JO*  IA4AJT100, Oil. **10*10*. 0910. 2. Oil. *92 

//  M>uivri*ticuu*c 

//ctiNne  woe 


// 

// 

// 

//• 

//• 


9ll9l*NUU9lkl.*IIIO*NUU9lll.VMlN*NUU9lll. 
2  NOW  *NUi  l  91 1 9 .  THtfl.O .  It  SHOT  *NUU9  til » 

6*  A9H*G0AJ0M.  191-f NC2C  IT 

9A9ANI2IAI  ON  TNI  IRIC  lIATf NINTl 


•00009010 
•000 00010 
•ooooooto 
00000000 
00000010 
900000*0 


//• 

//• 

//• 

NAN9 

02  9  AUl  2 

U1A69 

000 000  TO 

9R00 

A0AJI1 

9R06AAN  NAM 

000000*0 

//• 

11* 

9NM.9ROOI 

IIARAR2  HMfll  9I00IAN  9(91011 

OOOMIOO 

//• 

9119 

100 

NO.  09  CYLINOfAl  90*  92029001 

00000110 

//• 

OCR 

22*2 

910CRII1I  09  910)9  001.  9  2099  091,  ANO 

90000120 

//• 

920*9001 

OOOMDO 

//• 

99121 

NUU9I19 

OlNANl  90A  1112**2  2A9I  IN9U2 

00000140 

//• 

*9120 

NUL19III 

OSNANI  909  9IIIA91  2A9I  OU29U2 

OOOMIIO 

//• 

2MIN 

NU119KI 

OlNANl  909  YIM  HI  12092  IN9U2 

000001*0 

//• 

2NOU2 

Mill!  Ill 

OlNANl  909  YIM  HI  12092  OU29U2 

000M120 

//• 

2N19* 

0 

201*119  90*  IIM  HI  12092  IN9U2  2A9I 

OOOMIIO 

//• 

1VIU2I 

I22M9I 

OlNANl  909  92019001 

000001*0 

//• 

291 

ANY  YA9I  09122 

OOOMIOO 

//• 

001 

ANY  DUAL  Of Nil  22  YA9I  0*121 

OOOMIIO 

//• 

1219102 

Mill  *111 

OlNANl  909  9102  YA9I 

MO  90  220 

//• 

0MMI90 

— m»//cii 

II9C  9 CN* C9N06.il SION* 1 20AR .910RI 

0MM240 

//1T99CI0 

00  OI19*1h*.01N*11 19 

OOOMIIO 

UN  1 2a  21 14. 01N*C91Lf 9221 OUNI  .20t*K9*J09l0I.  OIS9* IH9.  •00000<»0 

OCMHI AACH2*  l.l  A9ll*l  . . . INI  00000220 

Of.  4.  I  *9C9N.99.lA9Cl*90.9kRII  21*0001  00090200 

UNI  2* 1 291 ..01919 1 .01 19*1  .CA216.0ILIII  I  .0IN*I9IS2Q.  *00000200 

OCA.HIAACHYM  00000900 

UNI  2*IVS0A,19ACf*ICYL.  1 19)1911  ,0IN*(IYIUY1,  *00000910 

QCflMC9M.y9t.l*[Cl*)»44.9lRlltl*C9lR,HIAACMV*ll,  *00000120 

OllfINU. 9*111  00000190 

UNI  fa  I  29$. .01999  1.0119*010  .0IN*t9ll2l ,  *000001*0 

UC9*NIARCMY*1  MO  00190 

nOHAMMlN  000001*0 

lYKlUf.A  00000120 

svsour**  oooMioo 

UNI 2*IY16A.19AC9*I Y9R.I 101 I <0IN*lt2tlNl •  *000001*0 

0CN*IMO^<H>.lRtCl*90.*lRlltl*)*00.HlAACMV-ll.  *00000400 

0119. IN)., 9*111  00000410 

01N*CS21INI.UN|l*tYiOA.t9ACf*<l9R,<IOtl.  *00000420 

0C4-INfC9N.9.0i.RSIII*l0.HIA9CHY*ll.0ll9*INlM.9Alll  00000410 

UNI  2a|21IIA.\9ACI*ir9R,  10.41  I  .OII9«INItl.Ofkl9lt.  *00000440 

OCN*HIA«C*<Y*t  00000490 

//C9I9102  IR9C  9IIN*Cf.*A9H.C0N0*l  1 4.LI.C11 1 .1 V1NI  000004*0 

//tirilll  00  0tl9.*H*,0IN*UI*l  00090420 

OI1*a  KILO.  0(1121 1  .OtN*9S2SUII.OC9*MIAACHY*l  00090400 

tVIOUIa*  000004*0 

UNI  2*  I 291,  ,02090 1 .01 19*010 .01N-CTH I N.  VOi*MR-(TH>IR.  *00000*00 

OCMMIAACHY*!  00000910 

UNIf *11*1.(029101. 01 S9«I(RI09. Of klTII.OINaRIHOUO,  *00000920 

OC9*|9IC9N*VNl«L9f Cl*)*44,0tRllll*lllR,H|A9CMY*ll  00000910 

OII9*im.O.OfllTfl.O|N*CIVIIN2,OeO*MIMCNV*l  00000940 

//9I0Y  2A99  00  llNI9*IQ09a.O99|9l.0tN*lt*l9l01,lA9fk*(,Ni),V0k*99IVATI,  *00000990 
//  0C4*HI  AACHY* 1  000009*0 

//  99N0  THIS  19  2N9  I  NO  09  INI  IN- 1 2*2*9  9MCI0UM  00000920 

//  9R9C  C919A0C 

//IN  00  • 

/* 


//9T0I900I  00 
// 

//• 

//9202900I  00 
// 

//9T0I900I  00 
// 

// 

//9TO49001  00 
// 

//9T01900I  00 
//f  20*200 1  00 
//9T0T900I  00 
//9TI0900I  00 
// 

// 

//92I 19001  00 
// 

//92I4900I  00 
// 


//9Y01900I  00 
//92049001  00 
//920I904I  00 
// 

//9204900I  00 
// 

//9TI0900I  00 


★REGION  also  shown  in  for  usa  with  HIARCHY  support.  If  OVERLAY 
la  used  or  entire  program  la  run  only  In  main  or  only  In  extended 
memory,  the  REGION  size  must  be  changed  accordingly. 


Figure  2-1.  Job  Control  Language  To  Run  AGAJ73  and  GDAJ07  Data. 
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3.  GENERAL  PROGRAMING  AIDS 


3.1  MARCO  PLOW  CHARTS 

The  flow  chart*  in  Figure*  3-1  through  3-6  describe  the  functional  struc¬ 
ture  of  the  program  without  regard  to  flow  by  subroutine.  Figure  3-1  shows 
the  total  program  structure  and  is  a  composite  of  Figures  1-1  through  1-13 
in  Volume  IX.  Figure  3-2  provides  some  detail  of  the' trim  process  in  Fig¬ 
ure  3-1.  Figure  3-3  amplifies  the  stability  analysis  shown  in  Figure  3-1. 
Figures  3-4  and  3-5  give  some  flow  logics  of  Runge-Kutta  Method  and  Ham¬ 
ming's  Predictor-Corrector  Method  employed,  respectively,  in  maneuver 
functions  in  Figure  3-1.  Figure  3-6  details  the  data  reduction  program. 

3.2  FORTRAN  SUBROUTINES  IN  AOAJ73 

The  FORTRAN  subroutines  contained  ity  AGAJ  are  listed  in  alphabetical 
order  including  the  main  program,  which  is ‘called  MAIN.  A  few  remarks  are 
made  for  each  subroutine  which  Indicate  its  general  purpose  or  use  in  the 
program.  In  the  case  of  subroutines  with  multiple  entry  points,  the  names 
of  all  the  entries  are  given  after  the  name  of  the  subroutine  in  the  order 
that  they  occur  in  the  subroutine. 

(1)  AJACOB.  AJACB1.  This  subroutine  handles  computation  of  quantities 
which  depend  upon  variables  which  are  changed  in  either  subroutine  TRIM 
or  subroutine  STAB  in  order  to  compute  partial  derivatives.  These  quanti¬ 
ties  are  then  calculated  and  used  in  the  computation  of  forces  and  moments. 

(2)  ALLMAT.  ALLVEC.  The  stability  analysis  uses  this  subroutine  to  com¬ 
pute  eigenvalues  and  "eigenvectors. 

(3)  ALSTAB.  This  subroutine  calls  ALLMAT  and  processes  and  prints  out 
the  results  of  the  stability  analysis. 

(4)  ANAL.  Output  of  this  subroutine  consists  of  the  total  summation  of 
forces  ani  moments. 

(5)  AZMI  NT.  It  Initialises  some  variables  which  will  be  used  by  sub¬ 
routine  AZMUTH. 

(6)  AZMUTH.  This  subroutine  in  the  rotor  analysis  does  the  calculation 
and  Integration  at  each  |>lade  azimuth. position. 

(7)  CDCL.  This  subroutine  uses  the  local  angle  of  attack  and  Mach  num¬ 
ber  plus  the  airfoil  aerodynamic  inputs  to  compute  the  steady-state  lift, 
drag,  and  pitching  moment  coefficients  for  each  blade  segment  of  the  rotor. 

(8)  CGZARM.  CCYARM.  CGXARM.  This  subroutine  calculates  moment  arms 
about  the  eg  for  the  aerodynamic  surfaces  whenever  eg  is  shifted. 


(9)  CLCD.  This  subroutine  is  similar  to  subroutine  CDCL  except  thet 
CLCD  computes  the  three  aerodynamic  coefficients  for  the  wing  and  stabili¬ 
sing  surfaces. 

(10)  CMC/LC,  This  subroutine  Interpolates  on  the  Carta  tables  to  pro¬ 
duce  the  contribution  of  unsteady  aerodynamics  to  the  pitching  moment.  It 
is  the  major  section  of  the  SUNS  unsteady  aerodynamic  model. 

(11)  CNTM.  AUXJ.  FLAP.  Some  of  the  forcing  functions  in  a  maneuver  may 
be  timed  to  start  after  the  rotors  have  been  stopped.-  This  subroutine  con¬ 
verts  those  relative  times  to  absolute  times.  It  also  varies  Jet  thrust 
and  RPM  dependent  flapping  stops  during  a  maneuver. 

(12)  CONSTB.  This  is  the  control  program  for  the  stability  analysis. 

(13)  CONTRM.  This  is  the  control  program  of  the  trim  segment. 

(14)  DAMPER.  RATI.  This  is  the  variable  damper  for  TRIM.  The  purpose 
of  this  is  to  gradually  dampen  out  oscillations  of  the  trim  Iterations. 

This  is  accomplished  by  checking  the  errors  generated  in  TRIM  against  an 
upper  limit  and,  whenever  all  errors  are  less  than  this  limit,  reducing 
both  the  partial  derivative  Increment  and  the  maximum  amount  which  one  of 
the  TRIM  variables  can  change  in  one  iteration.  This  second  entry  to  this 
subroutine,  RATI,  limits  and  applies  the  corrections  to  the  TRIM  variables. 

(15)  DATl.  This  block  data  contains  C. ,  C_,  and  CL.  tables  for  the  NACA 

0012  alTHTl.  L  0  M 

(16)  DAT3.  This  block  data  subroutine  contains  the  Carta  tables  used  in 
subroutine  CMCALC. 

(17)  DERIV.  This  subroutine  evaluates  the  derivatives  during  a  maneuver. 

(18)  DIFFER.  This  is  a  subroutine  of  numerical  differential  technique. 
Outputs  from  it  are  quantities  of  velocities  and  accelerations  required 
for  the  calculation  of  unsteady  aerodynamic  effects. 

(19)  ERRCHK.  This  subroutine  checks  possible  input  errors  in  the  program 
logic  group  iriiich  is  the  key  input  group. 

(20)  EXTORS.  It  recalculates  eg  location,  Inertias,  and  gross  weight 
when  any  external  store  is  dropped.  It  also  updates  aerodynamic  brake 
locations  if  brake  deployment  during  a  maneuver  is  exercised. 

(21)  FOCUS.  This  subroutine  is  in  the  rotor  analysis.  It  transforms 
fuselage  quantities  into  rotor  reference  and  calculates  accelerations  dur¬ 
ing  a  maneuver. 
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(22)  FUSFNM.  As  th«  turn*  implies,  this  subroutine  computes  fuselege 
aerodynamic  forces  and  moments.  It  slso  calculates  rotor  nacelle  drag 
contributions. 

(23)  FUSINT.  This  subroutine  converts  input  array  to  mnemonics  for  the 
fuselage  group.  It  also  calculates  eg  location  and  inertias  if  external 
stores  are  Included. 

(24)  OUST*  This  subroutine  is  entered  only  during  a  maneuver  in  which 
a  gust  is  bslng  genarated.  It  calculates  the  distance  of  each  part  of  the 
rotorcraft  from  the  start  of  the  gust  and  then  calculates  from  that  dis¬ 
tance  the  magnitude  of  the  gust  valoclty  at  each  point  on  the  ship. 

(25)  HAWS.  This  is  the  core  of  Hamming's  Prsdlctor-Corrector  Method. 

It  does  predictions,  checks  errors,  and  make#  corrections.  It  also 
handles  the  Integration  of  differential  aquations  and  calls  other  sub¬ 
routines  necessary  to  a  maneuver. 

(26)  HARM.  HARM1.  The  harmonic  analysis  for  blade  loads  at  the  trim 
point  is  done  by  this  subroutine. 

(27)  HPCG.  This  subroutine  calculates  the  first  three  time  points.  It 
is  used  along  with  trim  points  to  start  Hamming's  Predictor-Corrector 
cycles  in  subroutine  HAMMS. 

(28)  HRESP.  The  elastic  modes  during  the  quasi-static  trim  procedure 
are  handle;:  by  this  subroutine. 

(29)  IHBLD.  INBLDl,  INBLD2,  INBLD3.  This  subroutine  converts  input 
array  to  mnemonics  for  blade  related  data. 

(30)  INIT.  This  subroutine  initialises  the  arrays  for  the  printout 
during  a  maneuver  and  also  writes  those  arrays  on  disk  for  later  plotting. 

(31)  INRO.  The  function  of  this  subroutine  is  the  initialisation  and 
calculation  of  problem  constants  from  the  rotor  Inputs  only. 

(32)  INRTR.  This  subroutine  Initialises  some  of  the  rotor  related  data 
which  are  not  done  in  subroutine  INRO. 

(33)  INSCAS.  Initialisation  of  the  SCAS  Inputs  is  done  here. 

(34)  INSTAB.  This  subroutine  initialises  for  a  stability  analysis.  It 
calculates  the  partial  derivatives  needed  for  later  computing  the  frequency 
response. 

_  / 

(35)  INTFRQ.  It  interpolates  blade  natural  frequency  as  a  function  of 
rotor  rotational  speed  and  blade  collective  pitch. 
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(36)  INVERS .  Thi«  subroutine  calculates  tha  Inverse  of  the  mass  matrix 
before  the  call  to  ALLMAT. 

(37)  ICMAT.  This  subroutine  prints  tha  mass,  damping,  and  stiffness 
matrices  used  in  the  stability  analysis. 

(38)  ITRIM.  Included  in  this  subroutine  is  the  iteration  loop  of  the 
trim  section  of  the  program.  The  function  here  is  to  Iterate  to  a  trimmed 
flight  condition. 

(39)  ITROT.  This  subroutine  initialises  variables  for  subroutine  AZMUTH 
and,  when  specified  by  the  input  parameters,  controls  the  iteretlon  loops 
to  balance  the  rotor  flapping  moments.  It  also  calculates  the  trim  condi¬ 
tions  for  an  elastic  blade  . 

(40)  JACOBI .  As  the  name  of  this  subroutine  implies,  its  function  is  to 
calculate  the  Jacobian  for  use  in  the  Newton-Raphson  iteration  method  in 
TRIM  or  to  calculate  the  displacement  derivatives  for  use  in  the  stabil¬ 
ity  analysis. 

(41)  JFBCIN.  ATMINT.  This  subroutine  converts  input  arrays  to  mnemonics 
for  jets,  flight  constants,  bobwclght,  and  weapons  groups. 

(42)  JSTRED.  This  subroutine  reads  most  of  the  input  data  groupa. 

(43)  LGCINT.  The  program  logic  group  input  array  is  converted  to  mnemo¬ 
nics  in  this  subroutine. 

(44)  LIZE.  Initialisation  of  some  numerical  constants  is  done  in  this 
subroutine. 

(45)  LOADT.  This  subroutine  computes  and  prints  out  the  loads  on  an 
elastic  blade  at  the  trim  point. 

(46)  MAIN,  This  subprogram  reads  the  control  cards  which  direct  the 
flow  of  the  whole  problem.  The  path  is  selected  and  calls  are  initiated 
to  begin  working  the  problem.  Upon  return,  poaslble  errors  are  checked 
for  and  appropriate  action  is  taken.  If  an  error  is  detected,  an  error 
message  may  be  printed  out.  Then  the  program  either  terminates  execution 
or  starts  the  next  problem,  depending  on  the  severity  of  the  error. 

(47)  MANU.  This  subroutine  controls  the  time-variant  maneuver  segment 
if  Runge-Kutta  Method  is  used.  It  handles  the  Integration  of  the  differen 
tlal  equations  and  the  calling  of  the  other  subroutines  necessary  to  a 
maneuver. 

(48)  MANUV.  It  controls  the  time-variant  maneuver  segment  if  Hamming's 
Method  is  used. 
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(49)  MATRIX.  The  function  of  this  subroutine  Is  to  calculate  the  trans- 
formation  matrix  from  a  set  of  Input  Euler  ancles. 

(50)  MBAL.  This  subroutine  calculates  rotor  flapping  Increments  during 
each  iteration. 

(51)  MORPHS.  Damping  and  stiffness  matrices  for  a  stability  analysis 
are  calculated  here.  ; 

(52)  MNEM.  The  function  of  this  subroutine  is  to  calculate  problem  con¬ 
stants  from  input  data  and  to  initialize  for  a  problem. 

(53)  MODAL.  The  variables  which  are  functions  of  mode  shape,  frequency, 
and  mass  and  inertia  distributions  only  are  computed  in  this  subroutine. 

(54)  MODES.  This  subroutine  calculates  the  mass  matrix  for  a  stability 
analysis. 

(55)  MOMB.  SUBA.  SUBB.  This  subroutine  simulates  a  servomechanism  con- 
trolling  the  swashp late  while  the  main  rotor  Is  being  folded  horizontally. 

(56)  MTLT.  TTLT.  FLAT.  YAWP.  This  subroutine  handles  mast  tilt  during 

a  maneuver,  the  flat  tracking  mechanism,  and  the  yaw  control  servomechanism. 

(57)  NOPS.  N0PS1.  The  inputs  to  this  subroutine  are  the  number  of  azi¬ 
muth  stations  lor  use  in  the  rotor  analysis.  The  outputs  are  quantities 
which  are  functions  of  the  number  of  azimuth  stations. 

(58)  NPUTOT.  This  subroutine  prints  most  of  the  input  data. 

(59)  NSTED.  NSTEDl.  This  subroutine  calculates  some  variables  for  un¬ 
steady  aerodynamic  ef facts. 

(60)  PDZ1,  32.  Ths  inputs  of  this  subroutine  are  a  trim  partial  deriv¬ 
ative  matrix  (i.e.,  the  Jacobian)  and  the  typo  of  helicopter  or  rotor- 
craft  being  flown.  This  subroutine  then  changes  the  partial  derivative 
matrix  to  conditions  which  arc  known  to  hold.  Essentially,  this  is  an 
attempt  to  filter  numerical  "noise"  in  the  matrix. 

(61)  PUNCH.  It  punches  nonzero  elements  of  mass,  damping,  and  stiffness 
matrices  used  in  the  stability  analysis.  The  form  of  the  punched  output 
cards  is  explained  in  Volume  II. 

(62)  PYLINT.  It  converts  input  arrsys  to  mnemonics  for  the  dynamic 
pylon  group. 

(63)  QUAN.  This  subroutine  sets  mnemonics  from  the  integration  array  at 
the  end  of  each  time  point. 

(64)  RADBCN.  It  calculates  some  variables  used  by  subroutine  RADIAL. 
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(65)  RADIAL.  This  subroutine  In  the  rotor  analysis  does  the  calcula- 
tlons  end  integration  along  the  blade  radiua. 

(66)  RADOUT.  It  prints  output  from  subroutine  RADIAL  and  AZMUTH. 

(67)  RBDY.  It  initializes  radial  stationa  for  a  rigid  blade. 

(68)  READIN.  This  subroutine  contains  the  logic  for  reading  and  print¬ 
ing  the  input  data. 

(69)  REOATB.  It  handles  the  read-ln  of  airfoil  data  tablea. 

(70)  REDBMS.  It  handles  the  read-in  of  aeroelastic  blade  data  aa  well 
es  blade  mode  ahapc  data. 

(71)  REDCL.  It  reada  the  coafflcianta  of  lift,  drag,  and  pitching 
moment  o f  each  airfoil  data  table. 

(72)  REDID.  It  handles  the  raad-in  of  group  ID  cards. 

(73)  REDRWK.  RWKOUT.  It  performs  the  read-in  and  printout  of  rotor- 
induced  velocity  distribution  (RIVD)  tables. 

(74)  REDSWK.  SWKOUT.  It  performs  the  read-ln  and  printout  of  rotor 
wake  at  aerodynamic  surface  (RWAS)  tablea. 


(75)  RESTRT.  REST4,  RESTO.  REST 1.  REST 2.  Restart  tapes  are  written  or 
copied  by  this  subroutine.  For  the  first  case  of  each  run,  it  also  ini¬ 
tialises  all  variables  in  commons  to  sero. 


(76)  RCUST.  Subroutine  OUST  carries  wind  gust  calculations  to  the  rotor 
hubs.  This  subroutine  then  carries  the  calculations  to  the  blade  elements. 
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(77)  RMGKTA.  This  la  a  special  Runge-Kutta  routine  being  used  as  a 
starter  in  Hamming's  Predictor-Corrector  Method. 

(78)  ROTAN.  This  subroutine  may  bs  considered  to  be  the  outer  section 
of  the  rotor  analysis. 

(79)  RPTPO.  .It  controls  the  optional  trim  page  output. 

(80)  RTINIT.  This  is  the  control  routine  tdtlch  handles  the  initializa¬ 
tion  of  the  rotor. 

(81)  RTWAKE.  This  routine  calculates  blade  local  induced  velocity  when 
rotor  wake  table  option  is  used. 

(82)  SAVEIC.  RESTOR.  RE5TR1.  This  routine  saves  initial  conditions  dur¬ 
ing  calculation  of  the  first  three  time  points  of  Hamming's  Method. 


3-6 


jmj *>&  rtwmi&v* 


/ 


(83)  SAVOLD.  It  imi  the  values  of  some  variables  at  the  previous  time 
point  for  Hanning's  Method. 

(84)  S  CAS  IT.  The  highest  derivatives  in  the  differential  equations  for 
the  SCAS  are  calculated  here. 

(85)  SI  VAR.  This  subroutine  handles  the  initialisation  of  the  maneuver 
inputs  for  subroutine  VAR1  which  are  not  a  function  of  the  trio  point. 

(86)  SOLVE.  This  subroutine  solves  systems  of  linear  -equations  by 
Gaussian  elimination. 

(87)  STAB.  This  subroutine  computes  the  rate  derivatives  used  in  the 
stability  analysis. 

(88)  START.  The  function  of  this  subroutine  is  to  change  units  of  the 
input  arrays  and  set  them  equal  to  mnemonics  and  to  control  the  initializa¬ 
tion  segment. 

(89)  STBWAK.  This  subroutine  calculates  the  effect  of  rotor  wakes  on 
each  stabilizing  surfaces  when  a  surface  uses  RVAS  tables. 

(90)  STBZFM.  It  calculates  aerodynamic  forces  and  moments  at  all  stabi¬ 
lizing  surfaces. 

(91)  STBZIN.  The  function  of  this  routine  is  the  initialization  and  cal¬ 
culation  of  problem  constants  for  wing  and  stabilizing  surfaces. 

(92)  SUPERP.  This  subroutine  simulates  an  autopilot. 

(93)  SMAS.  This  subroutine  performs  the  function  of  linking  the  controls 
to  the  swashplates  with  the  appropriate  linkage  factors  and  phase  factors. 

(94)  SWSRAT.  It  calculates  some  intermediate  velocities  and  accelera¬ 
tions  . 

(95)  TABFIX.  This  subroutine  calculates  arrays  to  be  used  in  the  method 
of  calculated  entry  in  subroutine  TAB  I  NT, 

(96)  TABINT.  As  the  name  implies,  this  subroutine  does  a  table  inter¬ 
polation  for  C^,  0^,  and  tables. 

(97)  T ABOUT.  When  tables  are  used  for  C^,  CD>  and  C^,  this  subroutine 
prints  the  tables  after  the  input  data. 

(98)  TILT.  HASP.  TPFA.  This  subroutlna  controls  eg  shift  calculations 
for  several  different  manners  of  shifting  eg.  The  primary  function  is  in 
a  mast  tilt  maneuver.  It  provides  not  only  for  eg  shift  but  also  for 
changes  in  control  phasing  as  a  function  of  the  mast  tilt  angle.  Secon¬ 
dary  entries  handle  eg  shift  with  folding  of  a  rotor  either  when  it  is 
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being  folded  aft  after  being  tilted  forward  and  atopped  or  being  folded 
horicontally  after  a  atop. 

(99)  TIMLP.  Thie  aubroutine  containa  the  initialisation  neeaesary  at 
the  start  of each  time  step. 

(100)  TIVAR.  This  aubroutine  handles  tha  initialisation  of  the  maneuver 
Inputs  for  subroutine  VARI  which  are  a  function  of  the  trim  point. 

(101)  TRIM.  As  the  name  implies,  this  aubroutine  ie  the  primary  one  of 
the  section  of  the  program  for  finding  the  trimmed  flight  condition. 

(102)  TRMANU.  This  routine  sets  up  arrays  for  the  outputs  of  trim  as 
well  aa  maneuver  pages. 

(103)  TRMl.  The  primary  function  of  this  routine  is  to  initialise  vari¬ 
ables  used  by  time-variant  trim. 

(104)  TURN.  This  subroutine  handles  a  banked  turn.  Secondarily,  it 
handles  pushovers  or  pull-ups.  It  doss  so  by  checking  input  data,  pick¬ 
ing  up  proper  inputs,  and  doing  the  appropriate  initialisation  to  find  a 
trimmed  flight  condition. 

(105)  TVT.  This  aubroutine  controls  the  time-variant  trim  procedure. 

(106)  UNSTED.  This  is  the  major  section  of  the  UNSAN  unsteady  aero¬ 
dynamic  model  mentioned  in  Volumes  I  and  II. 

(107)  VARI.  This  subroutine  produces  the  effects  of  input  disturbancea 
during  a  tli-variant  maneuver.  The  Inputs  to  this  subroutine  are  the 
user-supplied  forcing  functions.  The  results  produced  by  these  functions 
are  the  output  of  this  subroutine. 

(108)  VCNTRL.  VCNT1.  VCNT2.  — .  VCNT10.  This  is  called  by  subroutine 
VARI.  It  calculates  the  effects  of  time-variant  maneuver  disturbance. 

(109)  VIND.  This  subroutine  calculates  the  Induced  velocity  of  a  rotor. 

(110)  WAG.  The  time-dependent  lift  change  by  the  Wagner  and  Kussner 
Method  is  computed  in  this  subroutine. 

(111)  WING.  As  the -name  Implies,  this  routine  computes  aerodynamic 
forces  and  moments  on  wings. 

(112)  WKTABN.  If  the  number  of  blade  radial  stations  input  to  rotor  wake 
table  is  not  twenty,  this  subroutine  Interpolates  the  table  into  twenty 
blade  radial  stations.  This  is  done  outside  the  Iteration  loops  so  that 

a  three-way  interpolation  can  be  reduced  to  two-way. 
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(113)  WRFM.  This  la  an  output  aubroutlna  which  wrltaa  out  tha  rotor 
force  and  nonant  summary  in  ahaft  rafaranca  and  tha  fuaalaga  rafaranca 
forca  and  nonant  aunmary. 


(114)  WRINST.  ICWRTM.  IOWRf 1.  Thla  aubroutlna  prlnta  output  during  tha 
conputation  of  partial  derivatives  for  a  atabillty  analyala. 


(115)  WRMANU.  Thla  aubroutlna  containa  tha  wrlta  atatananta  which  pro¬ 
duce  part  of  tha  manauvar  printout. 


(116)  WRMDAL.  WRMDL1.  WRMDL2.  Thla  aubroutlna  prlnta  tha  raaulta  of  eub- 
routlne  MODAL. 


(117)  WRMS.  Aa  a  raault  of  atabillty  analyala,  it  prlnta  out  mode  ahapaa 
aaaoclated  with  tha  rotorcraft  charactarlatic  roota. 


(118)  WROT1.  WROT.  Thla  la  another  output  aubroutlna  which  producaa  tha 
heading  for  the  printout  of  tha  Input  data  and  tha  trim  page. 


(119)  WRSTAB.  Thla  aubroutlna  prlnta  tha  rata  darivativaa  uaed  In  tha 
atabillty  analyala. 


(120)  WRTMNV.  Thla  aubroutlna  dafinaa  tha  output  arraya  for  trim  aa  wall 
aa  maneuver  pagea. 


(121)  WRTRIM.  Aa  the  name  lmpliea,  thla  routine  wrltaa  trim  page. 

(122)  WRVP.  Thla  la  atill  another  output  aubroutlne  which  producea  the 
printouta  of  the  partial  derivative  matrlcea  calculated  and  the  independent 
variablea  uaed  in  the  calculation  of  thoaa  derivativea. 

(123)  WSHDUF.  It  calculatea  fuaelage  effecta  on  downwaah  and  aldewaah 
anglea  at  winga  and  other  atablllxing  surfaces. 

(124)  XCONIN.  Initialization  of  all  control  linkages  Is  handled  by  this 
aubroutlne. 

(125)  XSTINT.  Thla  subroutine  converts  input  array  to  mnemonics  for 
external  stores/aerodynamic  brakes. 

(126)  XSTORE.  It.  calculates  aerodynamic  forces  and  momenta  at  each 
external  store/aerodynamic  brake. 

(127)  YRINIT.  This  subroutine  conditions  the  aerodynamic  inputs  for  the 
rotors. 

(128)  YSINIT.  This  aubroutlne  conditions  the  aerodynamic  inputs  for  the 
wing  and  stabilizing  surfaces. 


(129)  ZERO.  This  is  part  of  the  initialization  segment.  Every  variable 
in  this  routine  la  set  to  zero. 
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(130)  2LLCAL.  This  subroutine  computes  zero  lift  line  increments  at 
wings  and  other  stabilizing  surfaces. 

3.3  ASSEMBLY  LANGUAGE  SUBPROGRAMS  IN  AGAJ73 

AQAJ73  uses  four  assembly  language  subprograms  which  were  written  at  Bell 
Helicopter  Companyt  ABDUMP,  DATE,  DOTX,  and  PDSRED.  The  first  three  are 
available  from  the  SHARE  (an  IBM  user's  organisation)  FORTRAN  library. 
Although  PDSRED  itself  is  not  available  from  SHARE,  their  library  does 
include  two  routines  (OPEN  and  CLOSE)  which  are  generalisations  of  PDSRED 
and  can  be  adapted  to  replace  the  single  routine. 

(1)  ABDUMP.  This  diagnostic  routina  permits  the  program  to  terminate 
Itself  with  a  memory  dump  under  06/360.  It  does  nothing  except  issue  an 
ABEND  macro  with  uaer  code  0322.  The  routine  is  in  System/360  Assembler 
Language.  It  was  prepared  at  Bell  Helicopter  and  is  in  the  public  domain. 

(2)  DATE.  This  routine  with  argument  (NDATE)  returns  the  current 
system  date,  NDATE,  in  Gregorian  formt  mm/dd/yy.  NDATE  must  be  at  least 
eight  bytas  long.  The  routine  is  coded  in  System/360  Assembler  Language 
for  0S/36C  only.  It  uses  the  mecros  peculiar  to  the  MVT  option  of  0S/360. 
It  was  prepared  at  Bell  Helicopter  and  is  in  the  public  domain.  It  con¬ 
tains  the  following  entry  points t 

ENTRY  SETIME  (TINT).  This  entry  establishes  an  operating  time 
interval  against  which  to  check  program  operation.  This  interval  (TINT) 
is  in  minutes  in  floating  point  form.  The  routine  does  not  cause  execution 
to  terminate  at  the  end  of  the  designated  Interval.  This  entry  initialises 
TIMEX. 


ENTRY  TIMEX  (TU,  DT,  TL).  This  entry  checks  the  central  processor 
time  since  the  last  call  to  SETIME  or  TIMEX.  It  returns  three  argument 
values  in  floating  point  minutest 

TU  -  Time  since  initial  call  to  SETIME. 

DT  -  Time  since  last  call  to  TIMEX  or  SETIME. 

TL  -  Time  remaining  in  the  SETIME  interval. 

(3)  DOTX.  This  is  a  function  subprogram  with  the  following  argument 
list!  A,  IA,  B,  IB,  NP.  It  io  a  fast  method  of  computing  vector  inner 
product  using  double  precision  accumulation  of  sum  of  products  of  single 
precision  operands.  The  routine  is  physically  small,  about  twice  as  fast 
as  equivalent  FORTRAN  code,  and  quite  flexible.  It  provides  the  improved 
accuracy  of  double  precision  vector  products  without  the  increase  of 
storage  and  time  required  to  do  full  double  precision  analysis.  DOTX 
returns  a  double  precision  result  if  the  calling  program  types  it  as 
such.  The  arguments  are  described  as  follows  t 
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A  and  B  ara  tha  first  elementa  of  tha  two  alngla  praclalon 
vactora  to  ba  multiplied. 

IA  la  tha  nunbar  of  words  batwaan  auccaaalva  elements  of 

A.  It  Is  a  full  word  lntagar  graatar  than  aaro. 

IB  la  tha  number  of  words  batwaan  auccasslva  elements  of 

B.  It  Is  a  full  word  lntagar  of  any  valua. 

NP  la  the  number  of  element  products  to  ba  formed,  l.e., 

tha  numbar  of  elements  used.  It  la  a  full  word  Integer 
graatar  than  xaro. 

The  only  validity  check  made  on  tha  arguments  la  that  IA  and  NP  ba  graatar 
than  zero.  If  either  of  them  Is  out  of  range,  tha  routine  returns  a  zero 
result.  There  ara  no  other  error  exits. 


Tha  computations  of  the  subprogram  can  be  represented  by  the  following 
summation: 


X  *<NA  +  i*IA)*Z(NB  +  i*IB) 


(3-D 


where 


1-0 


X  la  the  output  of  the  subprogram, 
.th 


Y(NA)  la  the  NA  element  Ln  the  Y  array  (A  In  the 
calling  sequence), 

Z(NB)  la  the  NBth  element  ln  the  Z  array  (B  in  the 
calling  sequence),  and 

NP,  IA,'and  IB  are  as  defined  in  the  argument  list. 

For  example,  assume  that  tha  two  arrays  FORCE  and  DISP  are  dimensioned  to 
3  and  8  respectively.  Using  DOTX,  the  following  statement  could  be  pro* 
c cammed: 

WORK  -  DOTX(FORCE(l),2,DISP(2) ,3,3)  '  (3-2) 

The  equivalent  FORTRAN  statement  would  then  be  as  follows: 

WORK  -  FORCE(l)  *  DISP(2)  +  F0RCE(3)  •  DISP(5)  +  FORCE(S)*  DISP(8) 

The  routine  la  written  in  System  /360  Assembler  Language  for  any  level  of 
assembler.  It  was  prepared  at  Ball  Helicopter  Company  and  is  In  the 
public  domain. 

(4)  PDSRED.  This  routine  is  used  to  find  a  member  of  a  partitional 
dataset  and  to  make  It  accessible  to  a  FORTRAN  routine  through  normal 
sequential  READ  statements.  The  parameter  list  Is  as  follows: 
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ffltnbtr  an  8-byte  character  atrlng  which  idantifiad  tha 
daalrad  mambar 

&an  the  symbol  &  is  a  standard  punch.  Tha  symbol  sn  stands 
for  a  statement  number.  Zf  tha  mambar  is  not  found,  tha 
program  goes  to  this  ststement  in  tha  calling  routine 
rather  than  executing  tha  statement  following  tha  call 
to  PDSRED.  This  type  of  argument  is  an  IBM  FORTRAN 
language  extension. 

After  returning  from  PD5RED,  normal  FORTRAN  READ  on  unit  1  may  ba  issued 
and  an  and  of  dataset  is  sensed  as  with  any  other  dataset.  Before  another 
call  to  PDSRED  can  be  issued  for  the  same  PDS,  Re.  rind  1  must  be  issued. 

This  routine  is  coded  in  IBM  System/360  Assembler  Language  using  0S/360 
macro  facilities.  This  function  is  not  duplicated  in  any  other  operating 
system.  The  routine  is  compatible  with  all  0S/360  releases  through  and 
including  tha  21.6  release.  It  was  praparad  at  Bell  Helicopter  and  is 
in  the  public  domain. 

3.4  FORTRAN  SUBROUTINES  IN  GDAJ07 

The  FORTRAN  subroutines  contained  in  GDAJ07  are  listed  in  alphabetical 
order  including  the  main  program,  which  is  called  MAIN.  A  far  remarks 
are  made  for  each  subroutine  which  indicate  its  general  use  or  purpose 
in  the  program. 

(1)  CALCS 1.  This  subroutine  is  the  Interface  between  subroutine 
SCALIT  and  the  CALCOMP  plot  routines. 

(2)  CONPLT.  This  subroutine  controls  the  plotting  of  the  time  his¬ 
tories. 

(3)  CURVET.  This  subroutine  analyses  the  time  history  of  selected 
variables  during  a  maneuver.  This  analysis  la  accomplished  by  a  least- 
squares  curve  fit  followed  by  comparison  of  both  the  amplitude  and  phase 
angle  of  different  variables.  Then  one  variable  is  expressed  as  a  linear 
function  of  two  others. 

(4)  C81L.  The  function  of  tins  subroutine  is  the  transfer  to  a  disk 
of  maneuver  time,  history  data  which  has  been  stored  on  a  tape  or  the 
transfer  to  a  tape  of  maneuver  time  history  data  which  has  been  stored 
on  a  disk. 

(3)  DATP.  This  first  block  data  subroutine  contains  part  of  the  head¬ 
ings  for  plotted  time  histories. 

(6)  DAT1.  This  second  block  data  subroutine  contains  port  of  the  head¬ 
ings  for  plotted  time  histories. 


3-12 


f* 


(7)  FSPT.  This  subrout ins  controls  ths  harmonic  analysis  of  a  tims 
history. 

(8)  HARM.  This  is  the  harmonic  analyais  subroutine  used  by  subroutine 
FSFT. 


(9)  HEDING.  This  subroutine  generates  the  labels  for  the  time 
histories  from  the  data  stores  in  DATP  or  DATl. 

(10)  MAIN.  This  is  the  control  program  for  0ADJ07. 

(11)  PL0T1R.  This  subroutine  does  the  CALCOMP  plotting  of  the  results 
of  the  harmonic  analysis. 

(12)  PPLOT.  Ihls  is  the  printer  plot  routine  which  produces  plots  of 
tims  histories. 

(13)  SCALIT.  This  subroutine  sets  up  the  arrays  for  the  time  history 
plots. 

(14)  SCLFIX.  This  subroutine  calculates  scale  factors  for  the  time 
history  plots. 

(15)  WR0T1.  WROT.  This  subroutine  prints  the  headings  on  the  printer 
plots. 

3.5  LABELED  COMMONS  IN  AOAJ73 


There  arc  21  labeled  C0M40NS,  but  no  blank  COMMON,  in  AGAJ73.  Each  of  the 
C0M10NS  is  listed  below.  Any  special  ordering  of  variables  and  the  reason 
for  doing  so  are  given  along  with  some  general  comments. 

(1)  ANDO  IT.  The  first  9  variables,  HFRC  through  YSHRN,  in  this  COMON 
are  double  precision.  It  is  used  only  in  the  general  purpose  subroutines 
controlled  by  subroutine  ANAL. 

(2)  FLEX.  It  contains  most  of  the  variables  used  in  elastic  blade 
model  analysle. 

(3)  FORCMC.  This  COMMON  contsins  the  Carta  tables  used  by  subroutine 
CMCALC.  Since  it  is  defined  in  block  data  format  in  subroutine  DAT3,  it 
is  not  marked  as  a  COMMON-  block  in  Table  4*1. 

(4)  FORMIC.  This  COMMON  contains  most  of  the  variables  use<}'  ip  comput¬ 
ing  the  rotor-induced  velocity  distribution  from  the  table  stored  in 
FORWKl. 

(5)  FORWKl.  This  is  the  set  of  rotor-induced  velocity  distribution 
(RIVD)  tables  used  by  subroutine  RTVAKE. 
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(6)  FORY.  There  is  no  special  ordering  of  varlablaa  In  this  C0M40N. 

It  eonalata  of  tha  array  "Y"  oparatad  upon  by  Runge-Kutta  and  la  usad 

in  tha  Initialisation,  trim,  and  manauvar  aagmanta. 

(7)  FORYYY.  It  contains  most  of  tha  varlablaa  spaclflcally  usad  by 
Hanning's  Predictor-Corrector  Mathod. 

(8)  FORYYl.  This  COMMON  has  only  ona  array.  It  ratalns  all  tha  values 
of  previous  time  points  and  their  derivatives  used  by  Hanning's  mathod. 

(9)  FOSWK.  This  COMMON  contains  most  of  tha  variables  usad  in  com¬ 
puting  the  affects  of  the  rotor  wake  at  the  aerodynamic  surfaces  from 
the  tables  stored  In  FSWKl. 

(10)  FOSUKl.  This  is  the  set  of  rotor  wake  at  tha  aerodynamic  surfaces 
(RHAS)  tables  used  by  subroutine  STBWAK. 

(11)  INSTAR.  This  block  contains  most  of  tha  input.  It  is  used  in  the 
initialisation  segment. 

(12)  MAHAL.  Tha  first  59  variables,  XF  through  NQTR,  in  this  COMMON 
are  ordered  to  allow  I/O  and  other  manipulations  to  be  dona  on  an  equiva¬ 
lent  array.  The  next  11  variables,  COLSTK  through  BlTR,  are  ordered  for 
equlvalenclng  to  an  array.  Not  more  than  10  of  these  variables  are  uaad, 
and  the  array  KVAR  is  used  as  a  pointer  vector  to  choose  which  ones  end 
their  order.  The  next  22  variables,  AIM  through  PYLMD,  are  ordered  for 
equlvalenclng  to  the  array  VAR  in  subroutine  STAB  for  the  calculation  of 
derivatives.  The  variables  TAXL  and  TAXR  are  equiva  enced  to  an  array  fn 
subroutine  CNTM. 

(13)  STAMAN.  Variables  in  this  block  are  mostly  used  in  Initialization 
and  maneuver  segments.  Arrays  SHPGRP  through  SFTGRP  are  ordered  to  allow 
I/O  and  other  manipulations  to  be  done  on  an  equivalent  array. 

(14)  STAR  AD.  Most  of  the  variables  here  are  used  in  the  Initialisation 
and  general-purpose  segments. 

(15)  STARAN.  This  COMMON  is  used  in  the  initialization  and  general- 
purpose  segments. 

(16)  STBD.  This  COMON  is  used  only  in  the  stability  analysis. 

(17)  STRIAB.  COMMON  STRIAB  is  used  in  the  initialization,  trim,  and 
stability  analysis  segments. 

(18)  STRIMA.  The  first  16  variables,  TZM  through  TCLOCK,  are  ordered  for 
equlvalenclng  in  subroutine  MOMB.  The  next  six  arrays,  SGASRF  through 
SCASYC,  are  ordered  for  equlvalenclng  in  subroutine  INSCAS.  This  C0IM0N 

is  not  used  in  the  rotor  analysis. 
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(19)  TAB.  This  COMMON  contains  the  CQ  tabla  for  the  NACA  0012  airfoil. 

Since  it  is  defined  in  block  data  forostuin  Subroutine  DATl,  it  is  not 
marked  as  a  COMMON  block  in  Table  4-1. 

, 

(20)  TAB1.  This  COMMON  contains  Ct  and  Cn  tables  for  the  NACA  0012  air¬ 
foil.  Like  TAB,  it  is  block  data  in  DATl.  Consequently,  it  is  not 
identified  as  a  COMMON  block  in  Table  4-1. 

(21)  TOP  LOT.  This  COMMON  is  used  in  all  segments. 

3.6  LABELED  COMMONS  IN  ODAJQ7 

There  are  6  labeled  COMMONS,  but  no  blank  COMMONS,  in  GDAJ07.  Each  of  the 
COMMONS  is  listed  below,  together  with  pertinent  comments. 

(1)  INPLOT.  This  COMMON  is  used  by  subroutine  SCALIT  and  the  other 
subroutines  in  the  segment  for  plotting  time  histories. 

(2)  MAXMIN.  It  contains  the  maximum  and  minimum  values  of  the  speci¬ 
fied  variable.  It  is  primarily  usad  to  determine  the  scale  of  the  plot. 

(3)  PLOTD.  This  COMMON  contains  the  data  in  the  block  data  subroutine 
DATP  which  is  used  by  subroutine  HEDING  to  furnish  alphameric  headings 
for  time  histories. 

I 

(4)  PLOTDl.  It  has  the  data  in  the  second  block  data  subroutine  DATl 
which  is  used  by  subroutine  HEDING  to  supply  headings  for  time  histories. 

(3)  TIMPTS.  This  COMMON  is  chiefly  used  as  local  storage  in  sub¬ 
routines  CURVET  and  C81L. 

(6)  TOPLOT.  This  COMMON  contains  control  variables  and  is  not  the 
same  as  COMMON  TOPLOT  in  AGAJ73. 
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4.  DETAILED  PROCRAMMING  AIDS 
4.1  CONTROL  SECTION  CROSS-REFERENCE 


The  Control  Section  Cross-Reference  List  for  AGAJ73,  Table  4-1,  shows 
most  of  the  control  sections,  Including  C0M40NS,  which  are  referenced  by 
another  control  section.  The  exceptions  are  listed  on  the  first  page  of 
the  table.  These  exceptions  are  system  routines,  and  their  presence 
would  not  contribute  to  the  cross-reference  list.  The  remaining  pages 
of  Table  4-1  contain  the  control  sections  in  alphabetical  order  in  a 
column  on  the  left  side  of  the  page.  To  the  right  of  each  control  section 
name  is  the  cross-reference  information.  "LENCTK"  is  the  else  of  the  sub¬ 
routine  or  COMMON  in  hexadecimel  bytes.  "CALLED  BY"  gives  the  name  of 
each  control  section  referencing  the  control  section  whose  name  is  in  the 
column  on  the  left.  "IS  USED  BY"  gives  the  name  of  each  control  section 
referencing  a  control  section  in  the  "CALLED  BY"  list  or  by  another  con¬ 
trol  section  in  the  "IS  USED  BY"  list.  "CALLS"  gives  the  name  of  each 
control  section  referenced  by  the  control  section  whose  name  is  in  the 
column  on  the  left.  "USES"  gives  the  name  of  each  control  section 
referenced  by  a  control  section  in  the  "CALLS"  list  or  by  another  con¬ 
trol  section  in  the  "USES"  list. 


The  information  in  the  Control  Section  Cross-Reference  List  is  sufficient 
to  construct  the  sequence  of  subroutine  calls' from  which  an  overlay  struc¬ 
ture  can  be  made. 

As  noted  in  Section  3. 2, several  subroutines  have  multiple  entry  points. 
However,  the  Cross-Reference  List  (Table  4-1)  includes  only  the  primary 
names  of  subroutines;  it  does  not  include  the  names  of  any  of  these 
additional  entry  points.  In  the  case  where  a  call  to  a  subroutine  is 
actually  a  call  to  an  additional  entry  point,  the  primary  name  of  the 
subroutine  which  contains  the  specified  entry  point  is  used  in  the  Cross- 
Reference  List.  For  example,  Table  4-1  indicates  that  subroutine  AZMUTK 
cells  RADOUT  when  AZMUTH  actually  calls  AZMOUT  (a  second  entry  point  to 
RADOUT).  In  Section  3.2  the  names  of  any  additional  entry  points  follow 
the  name  of  the  subroutine;  however,  in  reading  the  program  listing,  it 
is  generally  more  useful  to  know  the  name  of  the  subroutine  which  con¬ 
tains  a  particular  entry  point  than  to  know  the  additional  entry  points 
of  a  subroutine.  Hence,  a  cross-reference  of  the  additional  entry  points 
to  the  subroutine  which  contains  them  is  given  in  Table  4-2. 


Table  4-3  contains  the  Control  Section  Cross-Reference  List  for  GDAJ07. 
It  is  read  and  used  in  exactly  the  same  manner  as  Table  4-1. 

4.2  PROGRAM  LISTINGS 

Listings  of  the  three  programs  supplied  under  this  cont?:act  are  on  file 
at  Eustls  Directorate,  USAAMRDL: 
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(1)  AOAJ73 

(2)  GDAJ07 

(3)  AS812A 

Not*  that  Program  AS812A  la  written  In  PL/1  rathar  than  FORTRAN. 

Program  AS812A  la  indepandant  of  th*  AGAJ73/GDAJ07  combination.  Th* 
purpoaa  of  the  program  la  to  parform  least-squared-errors  curve  flta  of 
wind  tunnel  data  which  provide  Inputa  for  the  Nominal  Angle  Equationa  of 
AGAJ73.  Th*  card  format  for  the  wind  tunnel  data  input  to  AS812A  waa 
choaen  to  be  compatible  with  th*  format  of  data  auppliad  by  the  LTV  Low 
Speed  Hind  Tunnel,  Delia* ,  Texaa.  Th*  data  la  read  in  by  the  atatement 
on  CARD  500  for  the  aequenced  aource  deck.  The  forr't  may  be  changed  by 
th*  progranmer  to  read  most  any  available  data. 

The  printed  output  of  AS812A  Uata  the  coefficients  of  the  equation  to 
which  the  data  has  been  fitted.  The  output  aleo  lists  the  Inputa  to 
AGAJ73  In  the  proper  sequence  and  format  with  the  ACAJ73  card  number. 
Punched  card  output,  for  direct  input  to  AGAJ73,  la  an  option  In  AS812A. 
Th*  job  control  cards  for  AS812A  are  shown  in  Figure  4-1. 

4.3  PROGRAMMING  OPTION  IN  NUMERICAL  INTEGRATION 

As  mentioned  in  Section  2  of  Volume  I,  there  are  two  methods  of  numerical 
Integration  techniques  in  AGAJ73.  The  first  one  is  the  Runge-Kutta 
method  which  is  programmed  as  the  default  option.  The  second  is  Hanning's 
method  iriiich  can  be  activated  only  by  setting  the  logical  variable  HMD 
to  TRUE.  This  variable  is  located  in  subroutine  LIZE  and  is  set  _to_ 

FALSE  in  the  program  delivered  under  this  contract. 

4.4  REDUCTION  OF  PROGRAM  MEMORY  REQUIREMENTS 

There  are  several  very  large  arrays  in  AGAJ73.  The  alee  of  some  can  be 
easily  reduced  without  any  significant  decrease  in  the  capabilities  of 
the  program.  The  arrays  which  can  be  easily  reduced  in  size  are  shown 
in  Table  4-4.  The  table  also  includes  the  name  of  the  COtMON  to  which 
each  belongs,  the  full  and  reduced  dimensions  of  each,  and  the  approxi¬ 
mate  savings  in  memory  requirements  for  reducing  each  array.  Up  to  330K 
of  storage  can  be  eliminated  with  minimal  effort. 

Referring  to  Table  4-3;  arrays  WKRTR,  WKSTB,  and  AUX  are  independent  of 
each  other.  Their  sizes  can  be  cut  individually  or  concurrently.  How¬ 
ever,  arrays  CURVED,  CURVEL,  and  CURVEM  are  dependent;  if  one  of  them  is 
reduced,  all  of  them  must  be  reduced  at  the  same  time.  If,  after  the 
program  is  link-edited  in  load-module  form,  the  programmer  wishes  to 
reduce  the  size  of  any  array  noted  in  Table  4-4,  it  is  necessary  to 
include  a  REPLACE  card,  before  the  INCLUDE  OLDLIB  or  SYSLMOD  card  in  the 
linkage  step.  This  REPLACE  card  must  th-  \  contain  the  names  of  all  arrays 
which  have  just  been  reduced  in  size. 
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In  Table  4-4  there  a."  two  columns  under  the  heedlnge  of  array  name  and 
dimension.  The  "Outride  RESTRT"  column  givea  the  name(s)  and  dimension(a) 
to  be  used  in  the  associated  COWON  block  when  it  is  included  in  any  sub¬ 
routine  other  than  RESTRT.  In  subroutine  RESTRT  the  information  listed 
under  "In  RESTRT"  must  be  used.  The  different  names  and  single  dimen¬ 
sions  are  in  RESTRT  simply  to  facilitate  the  taek  of  setting  the 
COMMON  blocks  to  zero.  Each  array  included  in  Table  4-4  is  discussed  in 
more  detail  below. 

VKRTR  contains  the  sets  of  coefficients  for  the  Rotor-Induced  Velocity 
Distribution  (RIVD)  tables  discussed  in  Sections  2.3.3. 1.3  and  2.3.3. 2 
of  Volume  II.  If  the  size  of  this  array  is  reduced  as  shown  in  Table 
4-4,  the  input  to  IPL(5)  must  be  equal  to  zero  and  in  subroutine  RESTRT 
FRK1(66000)  and  DO  34  I-l, 66000  must  be  changed  to  FRKl(l)  and 
DO  34  1-1,1,  respectively. 

WKSTB  contains  the  sets  of  coefficients  for  the  Roto.  Wake  at  Aerodynamic 
Surface  tables  discussed  in  Section  2.3.4  of  Volume  II.  If  '.he  size  of 
this  array  is  reduced  as  shown  in  Tabla  4-4,  the  input  to  IPL(37)  must 
be  equal  to  zero  and  the  inputs  to  XWG(29)  through  XWG(32),  XSTBi(29), 
and  XSTBi(32)  must  all  be  less  than  or  equal  to  100.  XSTBi  stands  for 

the  1^  stabilizing  surface.  In  subroutine  RESTRT,  FSK1(9000)  and 
DO  36  1-1,9000  must  be  changed  to  FSKl(l)  and  DO  36  1-1,1,  respectively. 

CURVED  is  the  array  where  drag  coefficients  of  the  airfoil  data  tables  are 
located.  If  the  size  of  this  array  is  reduced  as  shown  in  Table  4-4, 
the  input  to  IPL(2)  must  be  less  than  or  equal  to  2,  TABL(6350)  in  sub¬ 
routine  RESTRT  must  be  changed  to  TABL(3026),  and  TITLE(8,4)  and 
CURVED(  1100,4)  in  block  data  DATl  must  be  changed  to  TITLE(8,1)  and 
CURVED( 1 100 ,  1 ) ,  respectively. 

CURVEL  and  CURVEM  are  the  arrays  for  the  airfoil  lift  and  pitching  moment 
coefficients  respectively.  If  their  sizes  are  reduced  as  shown  in  Table 
4-4,  the  input  to  IPL(2)  must  be  less  than  or  equal  to  2,  TABL1(5375)  in 
subroutine  RESTRT  must  be  changed  to  TABL1(2150),  and  CURVEL(  1100,4) 
and  CURVEM(575,4)  in  block  data  DATl  must  be  changed  to  CURVEL(  1100, 1) 
and  CURVEM(575,1),  respectively. 

AUX  is  an  auxiliary  array  used  in  the  Hamming's  method.  If  this  array  is 
reduced  as  shown  in  Table  4-4,  the  default  option  of  numerical  inte¬ 
gration  in  maneuver  (the  Runge-Kutta  method)  must  be  used.  Also,  in 
subroutine  RESTRT,  FYY1(3440)  and  DO  104  1-1,3440  must  be  changed  to 
FYYl(l)  and  DO  104  1-1,1,'  respectively. 

There  are  other  large  arrays  which  can  be  reduced  if  storage  requirements 
are  still  excessive  after  taking  the  steps  discussed  above.  However,  the 
necessary  steps  for  reducing  them  are  more  complicated  and  in  some  cases 
mean  deletion  of  a  majo-  program  option,  and  in  most  cases  the  savings  is 
not  worth  the  effort.  Vie  programmer  should  consult  with  the  contractor 
before  attempting  any  storage  reductions  beyond  those  outlined  in  this 
report. 
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4.3  SWITCH  FOR  DIAGNOSTIC  DATA  FROM  STAB 

In  Section  3.2  of  Volume  II,  IPL(34)  la  defined  ee  a  awitch  for  obtaining 
diagnostic  data  during  the  stability  analysis  (STAB).  Since  the  data 
gensrated  by  this  switch  are  not  of  general  interest  to  the  user,  but 
can  be  useful  to  the  progratnaer,  the  function  of  IPL(34)  is  discussed 
in  this  Programmer's  Manual  rather  than  in  Volume  II.  The  function  of 
the  switch  is  described  below. 

In  STAB  there  are  up  to  22  independent  variables  which  may  be  incremented 
in  the  process  of  computing  the  stability  (partial)  derivatives.  The 
number  of  variables  actually  incremented  depends  on  the  number  of  degrees 
of  freedom  which  the  user  has  activated.  (See  IPL(30)  and  (32)  in 
Section  3.2  of  Volume  II).  In  each  STAB  case,  IPL(34)  can  be  used  to 
print  out  the  following  data  resulting  from  one  of  the  variables  being 
incremented! 


at  each  blade  station  and  each  azimuth  location  for 
each  rotor  (i.e.,  IPRINT  in  subroutine  RADIAL  does 
not  equal  zero,  which  calls  RADOUT) 

(2)  rotor  and  pylon  moment  data  (i.e.,  CONDI  in  subroutine 
MBAL  is  greater  than  1.3,  which  causes  printout) 

To  generate  this  output  for  a  particular  increment,  IPL(34)  is  set  to 
the  value  shown  in  Table  4-3.  Further  information  about  the  variables  in 
this  table  can  be  found  in  Section  4.11.2.1  of  Volume  II.  Note  that 
locking  out  a  degree  of  freedom  doea  not  change  the  correspondence  shown 
in  Table  4-5  between  IPL(34)  and  the  variables.  Also,  it  ia  only 
possible  to  obtain  this  extra  printout  for  one  variable  in  each  STAB 
case.  To  obtain  the  printout  for  more  than  one  variable,  the  case  must 
be  rerun  for  each  variable  of  Interest  with  IPL(34)  set  to  the  appropriate 
value  in  each  repeat  run. 


etc.) 


(1)  blade  element  aerodynamic  data  (or, 
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TABLE  4-3.  Continued 
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TABLE  4-4.  REDUCTION  OF  AGAJ73  S' 
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IE  (K.EQ.5)  RETURN  IF  (K.EQ.2)  RETURN 

K-5  K=2 

IF  (YRR(18,I).EQ.5.  .OR.  IF  (YRR(18,I).EQ.2.  .OR. 

YSTBZ(18,I).EQ.5.)  GO  TO  26  YSTBZ(18,I).EQ.2)  GO  TO  26 
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\  TABLE  4-5.  STAB  DIAGNOSTIC  SWITCH  IN  AGAJ73 

IPL(34) 

Variable 

1 

PUS.  u 

2 

PUS.  w 

3 

PUS .  Q 

4 

FUS.  V 

5 

PUS.  P 

b 

PUS.  R 

7 

H.R.  F/A  PYLON  RATE 

8 

M.R.  LAT  PYLON  RATE 

9 

T.R.  F/A  P\LON  RATE 

10 

T.R.  LAT  PYLON  RATE 

11 

M.R.  F/A  FLAP.  RATE 

12 

M.R.  LAT  FLAP.  RATE 

13 

T.R.  F/A  FLAP.  RATE 

14 

T.R.-  LAT  FLAP.  RATE 

15 

M.R.  F/A  PYLON  DISP 

16 

M.R.  LAT  PYLON  DISP 

17 

T.R.  F/A  PYLON  DISP 

18 

T.R.  LAT  PYLON  DISP 

19 

M.R.  F/A  FLAP.  DISP 

20 

M.R.  LAT  FLAP.  DISP 

21 

T.R.  F/A  FLAP.  DISP 

22 

T.R.  LAT  FLAT.  DISP 
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Figure  4-1.  Job  Control  Language  to  Run  AS812A. 


